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Cosmic strings are expected to form loops. These can act as seeds for accretion of dark matter, 
leading to the formation of ultracompact minihalos (UCMHs). We perform a detailed study of the 
accretion of dark matter onto cosmic string loops and compute the resulting mass distribution of 
UCMHs. We then apply observational limits on the present-day abundance of UCMHs to derive 
corresponding limits on the cosmic string tension CpL. The bounds are strongly dependent upon 
the assumed distribution of loop velocities and their impacts on UCMH formation. Under the 
assumption that a loop can move up to a thousand times its own radius and still form a UCMH, 
we find a limit of Gp < 1 x 10“^. We show, in opposition to previous results, that strong limits 
on the cosmic string tension are not obtainable from UCMHs when more stringent (and realistic) 
requirements are placed on loop velocities. 

PACS numbers: 98.80.Cq 


I. INTRODUCTION 

Cosmic strings (see e.g. [1-4] for reviews) are topolog¬ 
ical defects present in many theories of particle physics 
beyond the Standard Model. They are lines of confined 
energy density, analogous to defects such as vortex lines 
in condensed matter systems like superconductors and 
superfluids. In all particle theories that permit cosmic 
strings, a network of strings forms during a phase tran¬ 
sition in the very early Universe. Causality arguments 
[5, 6] show that this network persists to the present time. 
Using cosmological observations to hunt for gravitational 
effects of the energy trapped in cosmic strings is there¬ 
fore a powerful way to probe particle physics beyond the 
Standard Model [7]. 

Observable signatures of cosmic strings are typically 
proportional to the mass per unit length /i of the string 
[1-4], which is in turn related to the energy scale 77 
at which the strings form (/i ~ rf' [8]). Searching for 
cosmological signatures of strings thus probes particle 
physics in a “top down” manner, excluding higher en¬ 
ergy scales most easily. This makes searches for cosmic 
strings highly complementary to terrestrial accelerator 
experiments, which search for new physics via a “bottom 
up” strategy. 

As topological defects, by definition cosmic strings can¬ 
not have ends. They must either exist as part of a net¬ 
work of infinite strings, or as closed loops. Because they 
are relativistic, a segment of infinite string will typically 
have a translational velocity of the order of the speed of 
light. The network of infinite cosmic strings follows a 
“scaling solution” wherein the correlation length, which 
describes the mean curvature and separation of string 
segments, grows linearly with time. This causes the con¬ 
tribution of the network to the energy density of the Uni¬ 
verse to remain constant. The analytical arguments for 
the existence of this solution (see e.g. [4] for a review) 
have been confirmed by numerical simulations [9-15]. 


The scaling of the long string network is maintained by 
the formation of loops when segments intersect, remov¬ 
ing energy from the network. Studying the distribution 
of string loops numerically is much more demanding than 
following the network, because a much larger hierarchy 
of scales needs to be followed. However, current results 
indicate that loops also follow a scaling solution (see re¬ 
cent references in [9-15] for numerical evidence, and [16] 
for some more recent analytical work). 

String loops act as seeds for the growth of density per¬ 
turbations in the matter surrounding them. If they are 
present early enough, and persist for long enough, they 
can lead [17] to the formation of so-called ultracompact 
minihalos of dark matter (UCMHs; [18-26]). UCMHs are 
distinguished from regular dark matter (DM) halos by 
the epoch at which they undergo gravitational collapse. 
Regular halos do not transition to the nonlinear regime 
of structure formation until z < 30, whereas UCMHs 
collapse shortly after matter-radiation equality, in isola¬ 
tion. At this time, the background density field is still 
cold, smooth and essentially featureless. This means that 
UCMHs form by almost pure radial infall [22], giving 
them a far steeper central density profiles than regular 
cold dark matter halos: p ex [27-29] rather than 

pex r~^ [30]. 

If DM can self-annihilate, the rate goes as the square of 
the particle density. The steep density profile of UCMHs 
therefore makes them excellent candidates for indirect 
detection of DM [23, 31-35]. Searches for gamma- 
ray sources with the LAT instrument aboard the Fermi 
gamma-ray space telescope lead to the strongest limits 
on the cosmological density of UCMHs [24]. Gravita¬ 
tional leasing [22, 36, 37], neutrinos [38, 39], reionization 
[40, 41] and diffuse photon fluxes at various wavelengths 
[24, 42, 43] provide supporting limits. With these lim¬ 
its and a proper understanding of how to calculate the 
UCMH yield from a particular scenario, it becomes possi¬ 
ble to use UCMHs to place limits on the spectrum of pri- 
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mordial perturbations [24, 32, 38, 44], non-Gaussianities 
[45] and cosmic strings [17]. Here we provide an im¬ 
proved treatment of UCMH formation around cosmic 
string loops, and the resulting limits on the string tension 
Gfi. 

In Section II we review the important aspects of the 
cosmic string loop scaling solution, before following the 
accretion of DM by loops in Section III, and the sub¬ 
sequent formation of UCMHs. In Section IV we apply 
those calculations to determine the fraction of DM in 
loop-induced UCMHs, including the effects of loop ve¬ 
locities. We derive limits on the cosmic string tension in 
Section V, then summarize in Section VI. For the most 
part we use natural units, with c set to 1 . 


II. STRING LOOP SCALING 

An important ingredient in investigating the formation 
and evolution of dense structures seeded by cosmic string 
loops is the number density of loops per unit loop radius, 
n(i?, t). The quantity n{R, t)dR is the number density of 
loops in the cosmic string network at a time t with radii 
between R and dR. 

In this paper we use a one-scale model for the distribu¬ 
tion of string loops [46, 47], according to which all loops 
of initial radius i?i form at the same time Here i?i 

and are linearly related by a constant a\ 


The numerical value of a must be determined from sim¬ 
ulations; we adopt a = 0.05 [15]. String loops warp and 
twist as they evolve, so the probability that any given 
loop is exactly circular at any point in time is virtually 
nil. It is common to introduce a parameter /3 = 1/R that 
relates the radius of a loop to its length 1. Deviations 
from circularity can then be accounted for by allowing /3 
to differ slightly from 27r (although for our final limits we 
simply set (3 = 27r). 

The scaling solution implies that a constant number 
N of loops are formed per expansion time per Hubble 
volume, meaning that the number density of loops at the 
time of their formation is 


where z{t) is the cosmological redshift, and Zi{R) = 
2 :[ti(i?)] is the redshift at the time that loops of radius R 
were created. Making use of the fact that t oc (1 
during matter domination and t oc (1 -|- during ra¬ 
diation domination, the number density of loops is 

n{R,t) = Na^p-h-^R-^, (4) 

for loops formed during matter domination, and 

n(i?, t) = (5) 

for loops formed during radiation domination (as evalu¬ 
ated at some time t > t^q, where Gq is the time of equal 
matter and radiation). 

Emission of gravitational radiation by loops leads to 
a reduction of R with time, at an approximately con¬ 
stant rate proportional to G/i. The actual rate has some 
distribution over a population of loops, according to the 
loops’ individual geometries and the corresponding rate 
at which they can each emit gravitational radiation - 
but there is a typical value associated with the typical 
loop oscillation frequency. Combined with the fact that 
smaller loops are formed earlier, this means that there ex¬ 
ists a radius Rdec below which loops will have typically 
decayed by a given time t, 

i?dec(t) = (6) 

where 7 is another numerical constant determined from 
simulations [48]. We will use 7 = IOtt. 

From this point onwards we make the “fast decay ap¬ 
proximation”, where we assume that loops decay essen¬ 
tially instantaneously. In this approximation R(t) = i?i 
from the time of loop creation right up to decay, so R 
and i?i become essentially interchangeable in all the equa¬ 
tions that we have written so far. Under the fast decay 
approximation, the impact of gravitational radiation can 
be neglected for all radii R > Rdec- Below this value, 
due to the fact that there is always a tail in the distribu¬ 
tion of oscillation frequencies, and therefore some loops 
that shrink in radius at every rate less than the typical 
one, the physical number density of loops remaining will 
be proportional to R. The density per unit loop radius 
therefore becomes independent of R [1, 2], 


n(i?i, ti) = Ntr^ = Na'^p-‘^Rr^. (2) 


N is another constant that must be determined by sim¬ 
ulations; we take = 40 [15]. 

Neglecting, for the moment, slow decay of loops by 
emission of gravitational radiation, the physical radius 
i? of a string loop remains constant as the Universe 
expands. The number density redshifts, so that for 
t > ti (R) the number density of loops of radius R is 


n{R, t) 


z{t) + 1 1 ^ 
.Zi(i?)-bl_ 


n{R,ti ), 


( 3 ) 


n{R,t) = n{Rdec,t) for R < Rdec, (7) 

regardless of whether the loop was formed during matter 
or radiation domination. 


III. ACCRETION OF DARK MATTER BY 
STRING LOOPS 

Accretion of cold dark matter by a cosmic string 
loop leads to a spiky DM distribution [49, 50] (see also 
[27, 28]). If the resulting structure undergoes gravita¬ 
tional collapse sufficiently early, a UCMH will result. To 
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determine the total fraction of DM contained in loop- 
induced UCMHs at the present day, we must study the 
accretion of DM by string loops in some detail. 

There are several pitfalls to navigate in doing this. It 
is not valid to simply naively take all loops at tgq, and 
apply the growth factor from linear perturbation theory 
to the initial mass of each loop, as some loops may de¬ 
cay before such a period of accretion could become effec¬ 
tive. It would be equally incorrect to just eliminate loops 
that decay before a certain time (e.g. the present day, or 
teq), as such loops may have accreted enough mass be¬ 
fore their decay to have already created nonlinear struc¬ 
tures - possibly even during radiation domination. Even 
if such loops do not induce nonlinear structures by the 
time of their decay, the velocity perturbation that they 
induce will persist, and continue to grow in time after tgq 
even if the loop has already decayed. That growth may 
eventually lead to gravitational collapse in time to create 
a UCMH, even though the loop decayed much earlier. 
Thus, to study formation of UCMHs from string loops 
we must perform a careful study of the accretion of DM 
onto string loops both before and after 4 q, and before 
and after their time of decay 


td 




( 8 ) 


To compute the accretion of mass by loops of radius R 
produced at time we use the Zel’dovich approxima¬ 

tion [51]; specifically, the spherical collapse model. We 
consider fluctuations that are initially isothermal, where 
the initial fluctuation is exclusively determined by the 
cosmic string loop. The loop is the source of gravitational 
attraction, and over time will lead to an inhomogeneous 
DM distribution. 

We focus on a spherical shell of matter of physical 
radius r about the center of the loop. We are inter¬ 
ested in the time evolution of this radius, as the shell 
moves towards the center of the loop. It is convenient to 
parametrize this radius in the form 


r{x) = a{x)b{x)Q , (9) 


where a{x) is the scale factor, Q is the comoving coordi¬ 
nate associated with the spherical shell at the initial time 
ti, and b{x) measures the difference between the motion 
of the spherical region in the presence of the string loop 
compared to how it would evolve under simple cosmo¬ 
logical expansion. Here a convenient parameter for time 
is 


a{t) _ Zeq + 1 

a(feq) z{t) + 1 ■ 


( 10 ) 


In the case of cold dark matter we can neglect thermal 
motion of matter, and in this approximation each shell 
will be characterized by the relative mass fluctuation pa¬ 
rameter $(r) = (5M/M(r), where M(r) is the total mass 
within the shell of initial radius r, and 6M is the mass 
fluctuation due to the string loop (which is independent 
of r). 


As derived in [52], in the Zel’dovich approximation we 
have the following equation of motion for b{x): 


a;(a;-l-l) 


(fb 

dx"^ 


-I- 1 -I- -a: 


db 

dx 


l-b$ 


-5 = 0 . ( 11 ) 


where the information about which shell we are consid¬ 
ering is hidden in the value of d) for that shell. 

The solutions of this equation depend on whether we 
are in the radiation-dominated phase t < or in the 
matter-dominated phase t > teq, and whether the loop 
has decayed or not. In the radiation-dominated phase 
the DM is a sub-dominant component of matter, so on 
the sub-Hubble scales that we are considering when we 
study accretion by string loops, we expect only logarith¬ 
mic growth of the fluctuations after loop decay. In the 
matter-dominated phase, we expect linear growth in x. 
Hence, we must study the solutions of Eq. (11) in various 
cases. We will in fact find that the nonlinear mass grows 
linearly in x even in the radiation phase, as long as the 
seed loop is still present. 

The size at time x of a compact object formed via ac¬ 
cretion about the string loop is given by the radius of 
the shell which is “turning around” at the time x. The 
turnaround time txA for a fixed shell is the time when 
f = 0. At that time, the shell of matter disconnects from 
the Hubble flow to collapse, forming a virialized clump. 
Adopting the coordinates of Eq. (9), the turnaround con¬ 
dition becomes 

b + X— = 0 . (12) 

dx 

There is a critical turnaround time (with an associ¬ 
ated critical redshift, Zc) after which a collapsing overden¬ 
sity will not contain sufficiently pristine material to form 
by radial infall, and so cannot form a UCMH. The pre¬ 
cise time at which the radial infall approximation breaks 
down and a collapsing halo can no longer be said to form 
a UCMH is still rather uncertain [24, 36, 45]. Certainly a 
‘latest collapse redshift’ of Zc ^ 1000 is a conservative and 
very safe choice [22], but UCMH formation down to red- 
shifts as low as Zc ~ 0(100) is not inconceivable [24, 45]. 
Here we will assume Zc = 1000 for our final limits on the 
cosmic string tension G^,. Adopting a smaller yet still 
plausible redshift would lead to improved constraints, be¬ 
cause later collapse redshifts allow progressively smaller 
perturbations time to collapse and form UCMHs. 

We must verify that txA for the innermost shell of DM 
around a loop occurs before this critical collapse time, 
if the loop is to be said to have seeded a UCMH. The 
turnaround time is affected by the loop decay time; the 
earlier a loop evaporates, the longer it will take for the 
turnaround to occur. 


A. Formation and accretion before teq 

(Xi < 1, X < 1) 

In this section we study the accretion of DM by a 
string loop during radiation domination. The analysis 
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of this section is applicable to all loops formed before 
teq, but different parts of the analysis apply depending 
on whether the loop decays before or after tgq- We carry 
out the calculation in two different regimes. The first 
regime covers the period from loop formation (xi) up to 
either loop decay (xd) or equality (xgq = 1), whichever 
is earlier. This treatment is all that is required during 
radiation domination for loops that decay after equality, 
as in that case it is valid right up to tgq- For loops that 
decay before equality (xd < 1), the solution for the first 
regime must be matched on to the solution for the second 
at X = Xd. The second regime extends from loop decay 
to equality (x = 1). In the fast decay approximation, 
4) = 0 in this regime, as the loop is absent. 


1. Regime 1 (x < Xd) 

It is well known from the study of accretion of cold 
dark matter onto point seeds [49, 50] that the innermost 
shells are the first to decouple from the Hubble flow, turn 
around and collapse back onto the seed mass. To see 
whether any mass shell turns around by a particular time 
for a string loop with initial radius i?, we must hence 
focus on the shell of initial radius R. We first express the 
matter overdensity inside this shell in terms of cosmic 
string parameters. As long as the loop decay can be 
neglected, the value of the overdensity parameter 4)(i?) 
is constant. Starting from 




SM 

M{Ry 


(13) 


where 5M = Mioop = rPR is the mass of the loop and 


M{R) = -7ri?^pDM(ti) 


(14) 


is the total mass of DM contained within the region of 
initial radius R (the loop radius) at the time of loop for¬ 
mation, we find 


$(i?) = $(xi) = (15) 

Here and we have used the value of x at 

the time of loop formation to label the innermost shell 
instead of the loop radius R. Also, we have made use 
of the background Friedmann equation of motion (after 
rescaling the DM density to the time tgq and using the 
fact that the total density at that time is twice the matter 
density), and defined the constant k by 

’ (16) 


i.e. 


K = 


2^7rGpDM(teq)teq 

% 


which falls in the range 1/4 < n < 4/9. 


(17) 


During radiation domination x is small, so the second 
derivative term in Eq. (11) becomes negligible. With 
initial conditions 5(xi) = 1 for xi ^ 1, the approximate 
solution of (11) is 


5(x) 3 = 1 - ^x4>(i?), (18) 

for X ^ Xi- Inserting the result from Eq. (15) into this, 
we obtain 

b{xy = 1 — —. (19) 

Xi 


Inserting this result in the turnaround equation Eq. (12) 
yields the following expression for the turnaround time 


2 $(i?) 

= 2~^ f^KXi{Gfj,)~^ 


( 20 ) 


for the innermost shell. 

Consider now shells outside the innermost one, with 
initial physical radius r > R. For those shells the value 
of 4) is reduced to 

‘^’(^ = (7) m )- (21) 

The mass that has turned around by ‘time’ x is 
4:7r Q 

M(x) = —pDM(ti)rTA(x) , (22) 

where rTA(x) is the initial radius of the shell turning 
around at time x. From the first line of (20) it follows 
that 

1 

24>(rTA) ’ 

so that 

rTAixf = 2xi?^4>(i?). 

From (21) the nonlinear mass that has accreted around 
the loop at time x is 

M(x) = —M{R) = 2xMioop, (25) 

a^TA 

where xta still refers exclusively to the turnaround of 
the innermost shell. 

It is interesting to note the linear growth in x, the 
same growth obtained in linear perturbation theory after 
matter-radiation equality. It is also interesting to note 
that if the loop survives until x = 1, the nonlinear mass 
which has collapsed at that point is exactly twice the 
loop mass. 

Now we consider under which conditions the result¬ 
ing halo will actually collapse in this regime, as all ha¬ 
los that collapse already during radiation domination are 
sure to lead to UCMHs. Collapse of the innermost shell 


(23) 

(24) 
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in Regime I requires that xta < min(a;d, !)■ First of all, 
we see that only loops formed at ti < a~^jGfj,teq decay 
before equality, which corresponds to radii 

^ (26) 

^eq 

First we deal with the case where decay occurs before 
equality (xd < !)• Noting that during radiation domina¬ 
tion td/ti = (xd/xi)'^, from Eq. ( 8 ) we obtain 

Xd = (27) 

Comparing the expressions (20) and (27), we see 

that if decay occurs before equality, the condition for 
turnaround in Regime I is independent of the value of 
Xi (and therefore also R). For such loops, collapse in 
Regime I can only occur if 

Gm > 2-^13-^a^^K^fl. (28) 

For loops that decay after equality (xd > 1), the con¬ 

dition for collapse in Regime I instead becomes 

G/r > 2 - 3 /?- 5 / 2 a 3 / 2 ^^^ j ^ (29) 

For values of G/r larger than these critical values, 
nonlinear dark matter clumps will have formed before 
matter-radiation equality about cosmic string loops, and 
have therefore formed UCMHs. For R corresponding to 
decay before equality and values of G/r smaller than Eq. 
(28), we must also study how DM accretion continues 
between loop decay and equality to determine if UCMHs 
might still be created during radiation domination, but 
in Regime II instead of Regime I. 

2. Regime II (x > Xd) 

In cases where the loop decays before equality, we must 
continue to evaluate the evolution of the clump in the 
regime ccd < a; < I. Once again, the small x approxi¬ 
mation is valid and we can neglect the second derivative 
term in Eq. (11). The difference here with the calculation 
in Regime I is that we set $ = 0, as the loop is absent, 
and use b(x) from Regime I (Eq. 18) evaluated at x = Xd 
as the initial condition. This leads to solution 

2b^(xd) - 3 x + 3b^(xd)x + 3xd 
~ ~~~ 2 + 3xd 

= 1 - i3/2)xd^{R){2 + 3a;)(2 -t 3a;d)”^ (30) 

Inserting this into the turnaround condition (Eq. 12) 
gives the turnaround time 

~ 2$(R) 3xd<i>{R) ~ 2 

for the innermost shell. Making the correction | | 

to the final constant in order to match solutions exactly 


between Regimes I and II at Xd,^ we see that only for 
Gfj, larger than the critical value Eq. (28) do the sec¬ 
ond two terms give a positive correction to the result 
for Regime I (Eq. 20). This indicates (as expected) that 
there is always some growth in Regime II, but the cor¬ 
rection is small, so we can see that the additional growth 
is minimal. The critical relationship between Gfj, and 
R that results from demanding that collapse happens 
within Regime II (i.e. Eq. 31 < 1) is 

i?Aeq< 2^ ■3-^rhiGli) 

-2Va-3/2^i/2j-i^-i(G^)3/2 

-h 2^° • (32) 

Shells outside the innermost one also exhibit the re¬ 
duced growth. Following the same treatment as for 
Regime I, these will turn around at 

rTA(x)^ = 2R^d>(R) ^ + (33) 

3 -I- 2Xj 

leading to logarithmic growth of the nonlinear mass 

32; -L 2 

M{x) = 2 xdMioop 3^^ _|_ 2 • 

Here we can see both the linear growth of Regime I (the 
prefactor of 2 xd), and the nonlinear growth after Xd (the 
trailing correction). For loops that decay before equality, 
the UCMH mass at equality will therefore be 

M{x = 1) = (35) 

If Xd > 1 on the other hand, from Eq. (25) we see that 
M(x = 1) = 2Mioop. (36) 

We close this subsection with a warning: one might try 
to argue that if a nonlinear clump forms by x = Xd < 1 
then this clump would seed linear growth in the nonlin¬ 
ear mass during the period Xd < x < 1 in the same way 
that the loop itself seeds growth between Xi < x < Xd. 
This argument is incorrect because the outer shells, al¬ 
though sensitive to the presence of the seed loop, are not 
sensitive (by Birkhoff’s Theorem) to a re-distribution of 
mass within their own radii. This redistribution therefore 
leads to a nonlinear core but leaves a relative underden¬ 
sity between the core radius and the shell in question. 

B. Formation before teq, accretion after teq 
(Xi < 1, X > 1) 

After equality x > I, so the higher-order terms in x 
dominate the equation of motion (11). The approximate 


^ This is needed in order to account for the terms we neglected in 
Regime I when we took the approximation Xi <g 1. 
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form of the equation then becomes 


X 


2 


<Pb 

dx"^ 


3 db 
2^di 


1 

2 





= 0 . 


(37) 


Because we are only interested in 0{<1) negative cor¬ 
rections from unity to &, we can also take the term lin¬ 
ear in X and linearise it in b about 6=1. Doing this, 
and recasting in terms of the comoving displacement 
Ab{x) = 1 — b{x), the equation of motion becomes 

2x^A6"-b3a;A6'-(2$-b3)A6 = $, (38) 


where primes indicates derivatives with respect to x. 


1. Regime III (xd < 1) 


For loops that have already decayed before equality, 
<& = 0, so the equation of motion to solve is 

2x^Ab'' + SxAb' - 3A6 = 0, (39) 

which has the solution 

A6(x) 

0 

-f-a;[3A6(l) + 2A6'(l)]. (40) 

5 

The initial conditions are given by A6(l) and A6'(l), as 
calculated during radiation domination after loop decay, 
i.e. in Regime II according to (30): 


A6(l) = 1 - 6(1) 


= 1 - 


15 


1- yXd$(i?)(2 + 3a:d)-^ 


1/3 


A6'(l) = -XdHR){2 + 3xd)-^b{l)-^. 


(41) 

(42) 


We know that |a;d4’(i?)(2-|-3a:d)“^ < 1, so we can see 
that A6(l) A6'(l), allowing the A6'(l) terms in Eq. 

(40) to be neglected. Because we are interested in x > 1 
in this regime, we can also neglect the decaying solution 
proportional to leaving 

Ab{x) « ^xA6(l). (43) 

0 


Inserting this into the turnaround condition (Eq. 20) 
gives 


ccd^(R) ^ ^ 25 \ 

2 + 3xd 9 xta V 6xta 36x^^ J 


(44) 


Because xta > 1? this is approximately 

2 “h 3xd 
~ 3xd4>(R)' 


(45) 


Again, this does not perfectly match onto the solutions 
for Regimes I and II at x = 1 and Xd = 1, due to the 


various large-x approximations we have made along the 
way, but agreement can be forced with a simple 0(1) 
correction (at the expense of accuracy at very large x), 
giving 


a^TA = 


2 3xd 

10xd4>(R)' 


(46) 


In terms of the radius of a shell turning around at time 
X, this gives 


rTA(a:)^ = I0xXd'h(i?)i?^(2 -|- 3xd) 


-1 


(47) 


and a UCMH mass of 


M{x) 


\0xX(i]\^\qqp 

2 “h 3x,± 


(48) 


To have UCMH collapse occur in this regime demands 
xta < Xc, so that with Eq. (46) the condition on Gfi is 


G/t > 2 • 10 ^(2-I-3xd)^X(. ^/3 


This translates to 


R \ 


1/2 


<Ig, 


20 xc /3®/2 



(49) 


(50) 


which corresponds to positive R only when G/i > 2~^ ■ 
10 “^x“^/3“®a^7/^K^, the value below which UCMHs 
cannot form in the asymptotic limit i? —> 0. 


2. Regime IV (xd > 1, x < Xd) 

Eor loops that have not decayed by the time of equality, 
calculating the accretion of matter in the period between 
equality and loop decay requires solving the full inhomo¬ 
geneous ODE (Eq. 38). In this case the initial conditions 
at X = 1 are given by the calculation of accretion before 
equality and before decay, i.e. Eq. (19) from Regime I: 

\ 1/3 

4>(i?)j , (51) 

A6'(l) = ^6(1)-=^. (52) 

We see again that A6(l) A6'(l). 

The general solution to a linear inhomogeneous ODE 
is the sum of the general solution to the corresponding 
homogeneous equation, and any specific solution to the 
inhomogeneous one. Often, the simplest way to obtain 
the solution to an initial value problem like ours is to 
find the solution to the homogeneous equation that sat¬ 
isfies the initial conditions, and choose the specific solu¬ 
tion to be the one obtained with trivial initial conditions 
A6(l) = A6'(l) = 0. This way, the sum of the two solu¬ 
tions is guaranteed to satisfy the true initial conditions. 


Ab{l) = !-[!-- 
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The solution to Eq. (38) with A6(l) = A6'(l) = 0 is 


Abi(x) = 


$ 


3 + 2$ 


r--l/4 r 




/4(1 _ ^- 1 ) 


+a;'^/^(l + $-i) 



(53) 


with $ = (25 +16$)^/^. In general $ ^ 1, so $ ~ 5 and 

“'<*)» + (s-i) 

The solution to the homogeneous equation is approx¬ 
imately given by the solution from Regime III (Eq. 43). 
This is once more because $ is small, so the 2$ + 3 in 
Eq. (38) can be approximated to 3, reducing the homo¬ 
geneous form of Eq. (38) to (39). The decaying solution 
and A6'(l) terms can also be neglected once more, for 
the same reasons as in Regime III. We can now see that 
because the inhomogeneous solution (Eq. 54) is propor¬ 
tional to $, it is much smaller than the homogeneous 
solution. The homogeneous solution thus dominates en¬ 
tirely, and 


Ab{x) 


^xA 6 (I) 

o 



(55) 


The turnaround condition again results in a polyno¬ 
mial in xta similar to Eq. (44), and the lower-order terms 
in xta can be neglected because ccta > !> giving 


xta 


5 


(56) 


Using this to calculate the resulting mass of UCMHs that 
collapse in this regime gives 


M{x) 


3j^A7[oop 

5 


(57) 


Again, we can slightly correct this expression to properly 
match onto the Regime I solution at a: = I, resulting in 


M{x) = 2xMioop. (58) 

Using the correspondingly corrected expression for 
xtA) and demanding collapse before some final x (de¬ 
noted Xi) gives 

G/i > —• (59) 

2 Xi 

Here Xf = min(a::d, Xc)', if decay happens before the latest 
allowed redshift of collapse for UCMH formation then 
Xf = x^, otherwise xi = x^- 


allowed redshift of UCMH collapse. To solve for the ac¬ 
cretion history in this regime, we would take the solu¬ 
tion at the end of Regime IV as an initial condition, and 
evolve further from a;d using the homogeneous form of 
Eq. (38). However, the solution we obtained for Regime 
IV ended up being dominated by the homogeneous so¬ 
lution anyway, so all calculations for Regime IV apply 
directly to this regime. The only difference is that we 
must use Xf = Xc for determining the condition for col¬ 
lapse to occur in this regime, even though a:d < Xc- This 
gives 

—• (60) 
Z Xq 


C. Formation and accretion after teq 

{xi > 1, a: > 1) 

For loops formed after equality, we need to solve the 
equation of motion in the large x approximation, start¬ 
ing from initial conditions imposed at the time of loop 
formation Xi. 


1. Regime VI (x < Xd) 


Here we start from the initial conditions b(xi) = 1, 
b'(xi) = 0 as in Regime I (i.e. Ab(xi) = Ab'(xi) = 0), 
and solve the inhomogeneous equation as in Regime IV. 
The solution is 


Ab{x) 


$ 

3 + 2$ 



3 

5 



- 1 


(61) 


The resulting turnaround condition (again discarding 
smaller powers of xta/xi as xta ^ Xi) is 


Xta 




(62) 


which leads to the condition 


Gfi > 


5/3 ^a'^fxKxf 

8a:f — 20a;i 


(63) 


for collapse to occur in this regime. Here Xf is again 
the end-point of the regime, which will be the earlier of 
the decay time Xd and the latest allowed time of UCMH 
collapse Xc. From Eq. (62) we can also find the total 
mass of UCMHs collapsing in this regime. 


M(x) 



Moop- 


(64) 


3. Regime V (xd > 1, x > Xd) 


2. Regime VII (x > Xd) 


This regime deals with the period after a decay that 
occurs during matter domination, but before the latest 


After loop decay, we need to continue from Xd under 
the approximation that $ = 0. This means solving the 














same version of the equation of motion as in Regime III 
(Eq. 39), using A 6 (a:d) and Ab'{xd) from Regime VI (Eq. 
61) as initial conditions. This gives 


Abix)=U- 
5 V Xd 


- 3/2 


[A6(xd) - XdA6'(xd)] 


+ -— [3A6(xd) + 2xdA&'(xd)] • (65) 

OXd 

As X > Xd in this regime, we can drop the decaying solu¬ 
tion once more. We can also discard lower-order terms in 
Xd/xi after substituting in A6(xd) and A6'(xd) because 
a^d > a^i, giving 


Ab{x) 


3x$ 


Xd 


- 1 


5xd(3 -I- 2$) \ Xi 
With this solution, turnaround takes place at 


a^TA = 


X; ^ - X. ^ 


1 1 
^ 3 


( 66 ) 


(67) 


Demanding that Xta < Xc for collapse to occur in this 
regime gives 


G/i > 


15,0 f^Kxfxd 


24xcXd — 24xcXi — 40xiXd ’ 
with a corresponding UCMH mass of 


M(x) = 


2 , -1 - 1 . 2 
-x(Xi - Xj ) - - 


M, 


loop • 


Comparing this expression for the mass at x = Xd, 


M(xd) = 


2 Xd 16 

5 Xi 15 


M, 


loop 5 


( 68 ) 


(69) 


(70) 


to the corresponding expression at the end of Regime VI, 


M(xd) = 


2 Xd ^ 

5 Xi 


Ml 


loop 5 


(71) 


we see that there is excellent agreement. To enforce an 
exact match, we can simply correct Eq. ( 68 ) to 


G/i > 

and Eq. (69) to 
M(x) = 


5/3 ‘^a^f^Kxfxd 

8xcXd — 8xcXi — 12xiXd ’ 


2 / -1 -i\ 9 

5-(-, - -d ) - l5 


D. Summary 


Ml 


loop • 


(72) 


(73) 


Here we summarize the salient expressions from 
Regimes I-VII, arranging them into a set of five different 
scenarios for loop formation and decay, and a scheme for 
determining the present-day UCMH mass Mq and which 


(if any) regime UCMH collapse occurs in, for any given 
combination of G/t and R. We also illustrate the result¬ 
ing regions in Fig. 1. 

The expressions delineating the different regions are 
given in terms of the four critical times: Xd, xp Xc and 
Xto. The first two follow directly from R, 


Xi = 


Xd = 


_ti_ 

teq 


P R 


Oi t. 


eq 


/3 R 

iGfJ. teq 


(74) 

(75) 


where a = 1/2 during radiation domination and a = 2j3 
during matter domination. From these two expressions, 
we see that when G/t > a/y ~ 10 “^, decay happens 
essentially immediately and UCMHs cannot form. This 
constitutes a hard upper limit to the range of values of 

that we will consider. 

The second two critical times (Xc and Xto) are indepen¬ 
dent parameters that must be chosen for a given anal¬ 
ysis. The latest allowed time of UCMH collapse (xc) 
is discussed at the beginning of this Section; we adopt 
Zc = 1000, which leads to Xc = 3.12. The ‘turn off’ time 
Xto is the time during late matter domination at which 
linear accretion onto an already-formed UCMH ceases. 
As in previous papers [23, 24], we take this to correspond 
to the time at which structure formation has progressed 
far enough to catch most UCMHs up in bound struc¬ 
tures and prevent further accretion from the cosmological 
background: Zto 10, corresponding to Xto = 284. 

Scenario A (xj < 1 and Xd < 1): 

• Mo = 10xtoXdMioop(2 -I- 3xd)“^ (Eq. 48) 

• Gfi< Eq. (50) 

no UCMHs 

• Eq. (50) < G/i < Eq. (32) 

collapse in Regime III 

• Eq. (32) < G/i < Eq. (28) 

=4> collapse in Regime II 

• G/i > Eq. (28) 

=> collapse in Regime I 

Scenario B (xj < 1 and 1 < Xd < Xc): 

• Mo = 2xtoMioop (Eq. 58) 

• Gti< Eq. (60) 

no UCMHs 

• Eq. (60) < G/i < Eq. (59; X{ = Xd) 

=> collapse in Regime V 

• Eq. (59; X{ = Xd) < G/t < Eq. (29) 

=4> collapse in Regime IV 

• Gfi> Eq. (29) 

=4> collapse in Regime I 

Scenario C (x; < 1 and Xd > Xc): 
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• Mo = 2a;toMioop (Eq. 58) 

• Gfj, < Eq. (59; xt = Xc) 

no UCMHs 

• Eq. (59; Xf = Xc) < Gfi < Eq. (29) 

=4> collapse in Regime IV 

• Gm > Eq. (29) 

=> collapse in Regime I 

Scenario D (xi > 1 and Xd < Xc)' 

• Mo = [\xto{xl^ - Xj^) - Mioop (Eq. 73) 

• Gm < Eq. (72) 

no UCMHs 


• Eq. (72) < G/i < Eq. (63; Xi = a^d) 
=> collapse in Regime VII 

• Gfj, > Eq. (63; xt = Xd) 

=> collapse in Regime VI 


Scenario E (a;; > I and Xd > Xc)' 

• Mo = [^XtoX^^ — l] Mioop (Eq. 64) 

• Gfj, < Eq. (63; xt = Xc) 

no UCMHs 


Gfj, > Eq. (63; xt = Xc) 
collapse in Regime VI 


IV. ULTRACOMPACT MINIHALOS FROM 
COSMIC STRINGS 

A. Cosmological density of UCMHs 


The fraction /ucmh = ^^ucmh/^^dm of dark matter in 
UCMHs today is then simply 

r t?max 

/ucmh = pBm/ n{B!M)Md{GiJi,B!)dB!. (78) 

In contrast to the physical loop number density however, 
these two expressions need not account for the reduction 
in the number density of loops that have decayed by to, 
as the calculations of the previous section ensure that 
the impacts of loop decay are already accounted for in 
the expressions for Mq. Therefore, n{R,to) here should 
be understood to simply be Eq. (4) or (5) (depending 
on when the loops were formed), without the correction 
below Rdec given in Eq. (6). For loops created during 
radiation domination, Mq oc if ccd < 1 and Mq oc R 
if a:d > 1. For loops created during matter domination, 
Mq cx R^/^. The contribution to /ucmh from loops de¬ 
caying after equality is therefore dominated by those that 
decay close to equality. 

However, the available observational limits [24] on the 
number density of UCMHs apply to each UCMH mass 
independently, rather than to integrated mass ranges. 
These limits have been derived assuming a single UCMH 
mass, and therefore effectively apply to delta-function 
mass spectra. They can be easily converted to dif¬ 
ferential limits on /ucmh simply by dividing by the 
UCMH mass that they apply to, giving an upper limit 
on d/ucMn/dMo as a function of R. We show these lim¬ 
its in Fig. 2, along with the theoretical predictions, for 
different loop tensions and radii. 

The differential contribution to the fraction of DM in 
UCMHs of mass Mo{R,Gii) for a given combination of 
R and G/i is given by 


The actual physical number density of loops with ra¬ 
dius greater than or equal to R at time t is the integral 
of the number density per loop radius, over all loop radii 
greater than R, 

pRmayi 

nphys(i?,to) = / n{R',to)dR'. (76) 

Jr 

Here i?max refers to some maximum loop radius, e.g. the 
radius of loops created at equality if the number den¬ 
sity of loops with a;i < 1 is sought. Eq. (76) gives 
nphys{R,t) oc n{R,t)R for i?niax R- Referring to Eqs. 
(4) and (5), we see that nphys(R, t) cx R~^ for loops cre¬ 
ated during matter domination, and npdys{R,t) cx 
for loops born during the radiation era. In both cases the 
physical density is a steeply-falling function of the loop 
radius. 

The total present-day cosmological density of the dark 
matter in UCMHs, from cosmic strings of radius greater 
than or equal to 77, is likewise given by 


^^ucmh(G^, R) = 




n{R',to)Mo{Gfj.,R')dR'. 

(77) 


d/uCMH 

dMo 


d/uCMH 

dR 


/ 


dA^o 


= /Odm^(^> h)Md{G^i, R) 


dMg 




(79) 


where, as in Eqs. (77) and (78), n{R,to) is the number 
density not corrected for loop decay. G here is an order 
unity correction factor, ranging from 2/3 to 3 depending 
on the values of R and G/r, and essentially dictated by 
the inverse of the scaling of Mq with R: 


G{xi < 1, a:d < 1) 
G{xi < 1 , Xd > 1 ) 
G(xi > l,Xd < Xc) 

G(xi > l,Xd > Xc) 


2 -j- 3xd 

3 -I- 3xd ’ 

1 , 

6 xto(xr^ - - 9 

2xto(ccr^ - 2;/"^) - 9’ 
6 xto — 15xi 
2 xto — 15xi 


(80) 

(81) 

(82) 

(83) 


We see now that unlike the total cosmological density 
of UCMHs, the contribution to the differential UCMH 
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FIG. 1: Values of Gfi and loop radii R for which UCMHs can form. The different regions are marked according to which regime 
UCMH collapse occurs in, as defined in Section III. Diagonal lines show the borders between the areas in which loop decay 
takes place during radiation domination, during matter domination but before the latest permitted redshift of UCMH collapse 
(zc, which we set to 1000), and after UCMH collapse. The vertical line shows the border between the regions where loops form 
during either radiation or matter domination. The horizontal line is the largest value of Gfi for which the approximation of a 
constant decay rate (8), and our derived limits, hold. 


fraction has approximately the same dependence on R 
as the physical UCMH number density, and is always 
dominated by smaller loops. The strongest limits on G/i 
are therefore expected to come from the smallest values 
of R, as the observational limits fall off less steeply at 
low mass than d/ucMu/dMo oc (Fig- 2; [24]). 

B. Effects of cosmic string velocities 

Our calculations up to this point have assumed that all 
cosmic string loops are created at rest. Recent numerical 
simulations (e.g. [15]) show, however, that string loops 
are typically born with translational velocities that are 
a significant fraction of the speed of light. These large 
velocities are induced by the relativistic velocities of the 
long string segments from which the string loops are split 
off. After a string loop forms, its physical velocity is 


redshifted. 

Accretion onto a moving loop is not spherically sym¬ 
metric, and is less efficient than accretion onto a static 
loop. We take into account loop velocities in an approxi¬ 
mate manner, using the spherical growth formulae of the 
previous sections provided that the physical distance the 
loop has moved by the time it decays is smaller than some 
multiple k of the loop radius. If the loop moves further 
than this, we assume that no UCMH forms. 

For an initial physical velocity Vi^ the distance a loop 
of radius R has moved by the time it decays is 

Ar(/i)=afe)^'(^) (84) 

Evaluating this expression for loops that decay dur¬ 
ing radiation domination, we find that in order to have 
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FIG. 2: Some example comparisons between observational limits on the UCMH abundance, and the abundance predicted on 
the basis of cosmic strings. Here we show limits and predictions in terms of the differential fraction of dark matter in UCMHs 
as a function of Gfi, for four example loop radii. For these examples, we assume that a loop can move a distance of lOOOi? 
during its lifetime and still form a UCMH, and adopt a 100 GeV DM candidate that annihilates into bb pairs with cross-section 
3 X 10“^®cm®s"^ 


Ar{R) < KR, 


Vi < K 




(85) 


where K is a constant whose value we will vary. 

We take the distribution of initial velocities in each of 
three spatial directions to be Gaussian, as this provides 
a good fit to the results of numerical simulations [15, 17], 
leading to a speed distribution 


P{v) 


21 / 2^2 

7rl/2(u2)3/2 




( 86 ) 


The rate of UCMH formation becomes suppressed by a 
multiplicative factor S, given by the integral of ( 86 ) from 


zero to Vi'. 

^ = 3^1/2 („2)3/2- (8"^) 

Here we have Taylor-expanded the result about Vi = 0 
and kept only the first term; this is a very good approx¬ 
imation for 0{K) < 10^. In principle the distribution 
( 86 ) should also be cutoff at the speed of light, but for 
sufficiently small (u^) (< 0 . 1 ), the correction this induces 
is negligible. We will use the same suppression formula 
for all loops, even those that decay in the matter phase; 
the error this induces is also negligible. 

We adopt the value (z; 2)^/2 — 9 3 ^ assuming that the 
loop velocity distribution is similar to the long-string one. 
This quantity has been reported to be as low as 0.15 [17] 
or as high as 0.67 [15]. Given this uncertainty, and the 
difficulty in choosing an appropriate value of K, we show 



































12 


results for a range of possible velocity suppressions. Be¬ 
cause K and are degenerate in their effects on S, 

we simply vary K for the purposes of this illustration, us¬ 
ing K = Xc, 10 and 1000. For iF = 1, there is essentially 
no limit, as in this case the assumption of a constant loop 
decay rate (8) breaks down. We also show results with 
no suppression {S = 1). 

Multiplying the differential contribution (79) by the 
velocity suppression factor S and substituting in the loop 
distribution (5), we obtain 


d/uCMH _ ^IG'rrGCNa'^ 
dMo 


( 88 ) 


where X = ateq/{(3R) for loops formed in the radiation 
dominated era and X = 1 for loops formed during mat¬ 
ter domination. Inserting the values from (85) and (87) 
into (79), we obtain the following final formula for the 
differential fraction of mass in UCMHs induced by string 
loops formed during radiation domination: 


d/uCMH 

dMo 


16{2Tr)^/^a^K^ NC G 


x(zeq + l)-^/"ln-^(:^) . (89) 


This scales as K^, as R and as (G/i)^/^ (modulo 
logarithmic corrections in G/r). 

Note that our treatment of loop velocities differs in 
one major respect from that of [17]: their corresponding 
expression for (85) is missing a factor of y/G^, which 
arises from a missing factor of ^/t^■ Here we account 
for this extra factor, drastically weakening our limits in 
comparison to theirs. 


V. BOUNDS ON THE COSMIC STRING 
TENSION 

The fraction of DM in UCMHs of various masses is 
constrained [24] by the absence of dark matter self- 
annihilation signals observed by the FermiMAT tele¬ 
scope (see [53, 54]). In this subsection, we will apply 
these observational constraints on d/ucMu/dMo to de¬ 
rive limits on the cosmic string tension, using the theory 
of loop-induced UCMH formation that we developed in 
Section HI. 

Depending on the UCMH mass, the strongest lim¬ 
its may come from Galactic point sources, extragalactic 
point sources or the contribution of UCMHs to the Galac¬ 
tic diffuse gamma ray emission. The limits are summa¬ 
rized in Fig 2 (red curves), where d/ucMn/dMo varies in 
terms of Gfi, which translates into a variation in Mq for 
a hxed value of R. Based on Section HI, we also show 
the predicted differential fraction of DM in UGMHs as 
blue curves in Fig. 2. A constraint on G/i can thus be 
obtained for an associated R value. In Fig. 2, we show 
that for a few example radii, the limiting value of G/r is 


at the intersection of the two curves; values of Gp. for 
which the blue curve exceeds the red curve are observa- 
tionally ruled out. Note that the scaling with G/x agrees 
with what is expected from (89). 

Here we carry out the usual calculation of the UCMH 
core radius [24] for each combination of R and G/x, tak¬ 
ing the larger of the core radii implied by annihilation 
and angular momentum of dark matter. We further sup¬ 
plement this selection with a comparison to R itself, also 
demanding that the core radius must be at least as large 
as the loop. This, taken together with the fact that the 
loop radius is by definition far smaller than the actual 
halo turnaround radius, ensures that the non-sphericity 
of the loop does not interfere with the radial infall ap¬ 
proximation. 

We assume a dark matter mass of 100 GeV, a canon¬ 
ical thermal relic annihilation cross-section of 3 x 
10 “^® cm^ s“^ and 100% annihilation into bb quark pairs. 
These are relatively typical parameters for WIMP dark 
matter, and at about the limit of what is currently al¬ 
lowed by GMB and gamma-ray searches for dark matter 
annihilation [55-58] . It is however worth noting that the 
cross-section could be lower, the mass higher, or the final 
states less conducive to gamma-ray production than we 
have assumed. Although this would further weaken the 
resulting limits on G/x, in general UGMH limits on cos¬ 
mological scenarios are not strongly dependent upon the 
particle nature of dark matter (so long as it can actually 
annihilate) [24]. 

We plot our final constraints on G/x in Fig. 3, as a 
function of the loop radius R giving rise to the limit. It is 
worth remembering that because of the scaling solution, 
loops with all radii will be present - so the resulting limit 
on G/x is simply the strongest limit available at any R. 
The different curves we plot in this hgure are based on 
different assumed values of the velocity suppression factor 
S, assuming that the loop can move a different number 
of times its own radius before UGMH collapse is rendered 
impossible. 

The shapes of the limit curves can mostly be under¬ 
stood by considering how the observational limits vary 
over R (and correspondingly over Mq). At the high mass 
(high R) end, extragalactic sources lead to the tight¬ 
est constraints. At intermediate masses, galactic point 
sources provide the best limits. At the lowest values of R, 
it is the diffuse gamma ray emission that leads to the best 
constraints. The kinks in the curves below R ~ 10“® kpc 
correspond to the transition from galactic point source 
limits to extragalactic sources. Ignoring velocity sup¬ 
pression, we see that except at very small R, the shape 
of the limit curve tracks the theoretical minimum value 
of Gfj, at which cosmic strings can form UGMHs, rather 
than the observational limits on the UCMH abundance 
(compare to Fig. 1). 

Clearly the 5-factor selected can significantly impact 
the constraints on G/x. The limits improve as loops are al¬ 
lowed to travel further during their lifetimes, as a greater 
fraction of loops are assumed to be able to form UCMHs. 
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FIG. 3: Upper limits on Gjj. as a function of loop radius, for 
different velocity suppression factors. The limits have little 
dependence on the adopted DM model; we assume a 100 GeV 
DM candidate that annihilates into bb pairs with cross-section 
3 X 10“^®cm®s"^ 


It is only when loops are allowed to travel nearly 1000 
times their radii before UCMH collapse that the bound 
on Gfi becomes stronger than the accepted upper limit 
from CMB physics of < 1.7x 10"^ [59]. At K = 1000, 
the best constraint is G/r < 1 x 10“^ for an associated 
UCMH mass of ~ lO^M©. These results are clearly in¬ 
consistent with those of [17], who find limits of a similar 
order to this when assuming just K = 1. 

Neglecting to consider the 5-factor entirely, equivalent 
to cosmic string loops either having no translational ve¬ 
locity or their velocities being irrelevant to their ability 
to form UCMHs, improves the upper limit on G^ by 
several orders of magnitude. The best constraint, at a 
UCMH mass of 3 x IQ-^Mq, is G/r < 2 x IQ-^^. Of 
course, the most physically rigorous treatment demands 
a velocity correction factor, so this case is not realistic. 
The results show how small differences in the assump¬ 
tions with regard to this aspect can propagate into order 
of magnitude differences in the constraint on G/r. Our 
results can therefore be considered to present limits rang¬ 
ing from “best case” to “worst case” scenarios. Improved 
knowledge of the loop velocity distribution and simula¬ 
tion of UCMH formation around a moving seed would 
significantly increase the precision of the bounds. 

The final bound on G/r is also, unfortunately, sensitive 
to the parameters that describe the cosmic string loop 
distribution, in particular the values of a and For the 
values we have used in this paper, drawn from [15], our 
best limit on G/i with e.g. K = 1000 is G/j, < 1 x 10“^. 
However, for a = 0.5, we would obtain a more strin¬ 
gent bound: G/r < 10“®. However, one cannot simply 


extrapolate the constraint obtained when increasing a 
without taking into account the corresponding increase 
in the lower value of G/r below which no UCMHs form 
(see e.g. 50) . Similarly, if we reduce /3 by a factor of ten, 
we also obtain a tighter bound: G/i < 2 x 10“^°. 

The largest uncertainty to our results, however, comes 
from the treatment of the effect of loop velocities. We 
have modelled the effects with a factor of K which we 
varied, K = Xc being the most restrictive assumption 
and larger values of K being less restrictive. For this 
issue, an improved analysis might be possible. For exam¬ 
ple, we could adopt a “delayed start approximation” in 
which instead of removing the effect of loops which have 
too large initial velocities, we follow the loops and let the 
accretion start at a later time when the velocity has red- 
shifted to a sufficiently low value (we thank the Referee 
for making this suggestion). Another way to improve on 
our present analysis would be to study the cylindrical 
accretion of matter onto a loop which is initially rapidly 
moving and to use the part of that mass which is suffi¬ 
ciently spherically distributed to yield an ultra-compact 
object. We leave this issue for followup work. 


VI. CONCLUSIONS 

Because many particle physics models beyond the 
Standard Model give rise to cosmic strings, it is inter¬ 
esting to explore bounds on the cosmic string tension /i 
from a variety of cosmological observations. In this paper 
we have considered what limits can be set on the basis of 
the non-observation of gamma ray signals from DM an¬ 
nihilation events in ultra-compact mini-halos (UCMHs). 

UCMHs would form by the accretion of cold dark mat¬ 
ter by cosmic string loops at high redshift. The number of 
such UCMHs increases as G/r decreases. Assuming that 
the DM is self-annihilating, we obtain bounds on the cos¬ 
mic string tension. These bounds depend very sensitively 
on the parameters describing the cosmic string loop dis¬ 
tribution, in particular the value of a, as well as on the 
parameters that describe the loop velocity distribution. 
Assuming that loops are still able to form UCMHs even 
if they move quickly enough to travel a thousand times 
their radii before the nascent UCMH collapses, we derive 
a limit of G/i < 1 x 10“^. Tightening this assumption 
however leads to much weaker limits. 

Our basic method is similar to that of [17], but our 
final results are not consistent with theirs, mainly due 
to a more careful treatment of the velocity suppression 
factor on our part. 

Here we have also developed the analytical theory of 
the accretion process in great detail, which we hope will 
prove a useful reference for future work on loop-induced 
UCMHs. Although we have shown that this method only 
yields improved constraints over the currently accepted 
upper limit on G/r when unrealistic assumptions are 
made about loop velocity distributions, we have signif¬ 
icantly improved the theoretical understanding of mini- 
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halo formation from cosmic strings, and the accuracy of 
gamma-ray limits on G/r. Our results also suggest that a 
better understanding of the velocity distribution of loops 
and the formation of UCMHs around moving seeds might 
be fruitful for helping to constrain G^. 

All of the numerical UCMH routines used in this pa¬ 
per, as well as those from earlier papers [23, 24, 45], can 
now be found in v5.1.2 of the public software package 
DarkSUSY (www.darksusy.org [60]). 
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